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Abstract 

Transition disk objects are pre-main-sequence stars with little or no near-IR 
excess and significant far-IR excess, implying inner opacity holes in their disks. 
Here we present a multifrequency study of transition disk candidates located in 
Lupus I, III, IV, V, VI, Corona Australis, and Scorpius. Complementing the in- 
formation provided by Spitzer with Adaptive Optics (AO) imaging (NaCo, VLT), 
submillimeter photometry (APEX), and echelle spectroscopy (Magellan, Du Pont 
Telescopes), we estimate the multiplicity, disk mass, and accretion rate for each 
object in our sample in order to identify the mechanism potentially responsible 
for its inner hole. We find that our transition disks show a rich diversity in their 
SED morphology, have disk masses ranging from < 1 to 10 Mjup and accretion 
rates ranging from < 10 -11 to 10~ 7 ' 7 M yr _1 . Of the 17 bona fide transition disks 
in our sample, 3, 9, 3, and 2 objects are consistent with giant planet formation, 
grain growth, photoevaporation, and debris disks, respectively. Two disks could 
be circumbinary, which offers tidal truncation as an alternative origin of the inner 
hole. We find the same heterogeneity of the transition disk population in Lupus 
III, IV, and Corona Australis as in our previous analysis of transition disks in 
Ophiuchus while all transition disk candidates selected in Lupus V, VI turned 
out to be contaminating background AGB stars. All transition disks classified as 
photoevaporating disks have small disk masses, which indicates that photoevap- 
oration must be less efficient than predicted by most recent models. The three 
systems that are excellent candidates for harboring giant planets potentially rep- 
resent invaluable laboratories to study planet formation with the Atacama Large 
Millimeter/Submillimeter Array. 

Subject headings: accretion disk -- binaries: general -- line: identification - 



planetary systems: protoplanetary disks — stars: pre-main sequence 



Introduction 



Low-mass pre-main-sequence (PMS) stars are generally separated in two different 
classes, accreting classical T Tauri stars (CTTSs) with broad Ha emission lines, blue 
continuum and near infrared excess and non-accreting weak-line T Tauri Stars (WTTSs) 



with narrow symmetric Ha emission lines (e.g. Bertout 1984). While CTTSs typically 
show large excess emission from the near-infrared to the millimeter, WTTSs often have no 
infrared (IR) excess at all. Only a relatively small fraction of T Tauri stars are observed 
in an intermediate transition state with little or no near-IR excess and significant far-IR 
excess. This clearly indicates that once the inner disk starts to dissipate, the entire disk 
disappears very rapidly ( Wolk fe Walter|[l996 Andrews fc Williams||2005 Cieza et aT||2007 ). 
The missing near-IR excess combined with the clear presence of an outer disk is the defining 
characteristic of transition disks. However, a precise and generally accepted definition 
of what constitutes a transition disk object does not yet exist. The most conservative 
definition of transition disks, often labeled classical transition disks, consists of objects with 
no detectable near-IR excess, steeply rising slopes in the mid-IR, and large far-IR excesses 
(e.g. |Muzerolle et~aTl|2006| |Sicilia-Aguilar et al.|[2006| |Muzerolle et~al]|2010[ ). Being less 



restrictive, objects with small, but still detectable, near-IR excesses (e.g. Brown et al. 2007 



Merin et al. 2010) can be included, until considering objects with decrement relative to the 



Taurus median Spectral Energy Distribution (SED) at any or all wavelengths (e.g. Najita 



et al. 2007 Cieza et al. 2010). Throughout this paper we follow the latter and broader 



definition. However, one has to be aware that this broad definition still is mostly sensitive 
to inner opacity holes but may overlook pre-transitional disks with a gap separating an 
optically thick inner disk from an optically thick outer disk. Such systems have been 
identified from Spitzer IRS spectra ( Espaillat et al.||2~0~07 ), but can be missed by photometric 
selection alone. 



The Spitzer Space Telescope generated a huge database containing IR observations of 
PMS stars in star-forming regions. Most importantly, Spitzer products such as the catalogs 
of the Cores to Disks (c2dj£j and Gould Belt Spitzer (GB) Legacy projects (Spezzi et al. 



2011 Peterson et al.|[2~011[ ) provide SEDs from 3.6 to 24 /xm for large numbers of PMS stars. 
One of the most interesting results concerning transition disk studies with Spitzer has been 



the great diversity of SED morphologies (see Williams & Cieza 2011, for a review). The 



widespread of IR SED morphologies found in transition disk objects cannot be adapted 
to the classical taxonomy to describe young stellar objects (YSOs) such as the Class I, 



II, III definitions from Lada (1987). Cieza et al. (2007) quantified the richness of SED 



morphologies in terms of two parameters based on the SED shapes considering the longest 
wavelength at which the observed flux is dominated by the stellar photosphere, A tU moff, and 
the slope of the infrared excess, a excess , computed from A tU moff to 24 /xm. 

Studying the diverse population of transition disks is key for understanding 
circumstellar disk evolution as much of the diversity of their SED morphologies is likely 
to arise from different physical processes dominating the disk's evolution. Evolutionary 



processes that may play an important role include viscous accretion (Hartmann et al. 



1998), photoevaporation (Alexander et al. 2006), the magneto-rotational instability (MRI) 



(Chiang <__ Murray-Clay 2007), grain growth and dust settling (Dominik <__ Dullemond 



2008), planet formation (Lissauer 1993 Boss 2000) and dynamical interactions between the 



disk and stellar or substellar companions (Artymowicz & Lubow 1994) 



As discussed by Najita et al. (2007); Cieza (2008); Alexander (2008), one can 



distinguish between some of these processes if certain observational constraints, in addition 
to the SEDs, are available. To this end, we are performing an extensive ground-based 
observing program to obtain estimates for the disk masses (from submillimeter photometry), 



1 http://irsa.ipac. caltech.edu/data/SPITZER/C2D/doc/c2d_del_document.pdf 
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accretion rates (from the velocity profiles of the Ha line), and multiplicity information 
(from AO observations) of Spitzer-selected disks in several nearby star-forming regions. 



Our recently completed study of Ophiuchus objects (Cieza et al. 2010, hereafter Paper I) 
confirms that transition disks are indeed a very heterogeneous group of objects with a 
wide range of SED morphologies, disk masses (< 0.5 to 40 Mjup) and accretion rates 
(< 10~ n to 10~ 7 M yr^ 1 ). Since the properties of the transition disks in our sample 
point towards different processes driving the evolution of each disk, we have been able to 
identify strong candidates for the following disk categories: (giant) planet-forming disks, 
circumbinary disks, grain-growth dominated disks, photoevaporating disks, and debris 
disks. 

We here follow the same approach as in Paper I in performing multiwavelength 
observations to derive estimates on disk masses, accretion rates, and multiplicity. We 
present submillimeter wavelength photometry (from APEX), high-resolution optical 
spectroscopy (from the Clay, and Du Pont telescopes), and Adaptive Optics near-IR 
imaging (from the VLT) for Spitzer-selected transition circumstellar disks located in the 
following star forming regions: a) Lupus: I, III, IV, V, VI, b) Corona Australis (CrA), and 
c) Scorpius (Sep). 

2. Transition disks in Southern star-forming regions 

The Lupus clouds constitute one of the main southern nearby low-mass star-forming 
regions containing the following sub-clouds at slightly different distances: Lupus I, IV, 



V, VI at 150 ± 20 pc and Lupus III at 200 ± 20 pc flComer6n||2008| ). The clouds are 
situated in the Lupus-Scorpius-Centaurus OB association spanning over 20 deg in the sky. 
Their population is dominated by mid M-type PMS stars, but some very late M stars or 



substellar objects have been found as well thanks to Spitzer capabilities (see Comer6n|[2008, 



for a review). In general, the ages of the Lupus clouds are estimated to be ~ 1.5 — 4 



Myr, (Hughes et al. 1994 Comeron et al. 2003). However, a comprehensive analysis using 



Spitzer IRAC and MIPS observations in combination with near-IR (2MASS) data has been 
performed for Lupus I, III, and IV by the c2d Legacy Project ( Merfn et al.||2008 ) and Lupus 



V, VI by the Gould Belt Legacy Project (Spezzi et al. 2011) and revealed a significant 



difference between the sub-clouds. While Lupus I, III, and IV are dominated by Class 
II YSOs, Lupus V, VI mostly contain Class III objects. This has been interpreted as a 



consequence of Lupus V, VI being a few Myrs older than Lupus I, III, and IV by Spezzi 



et al. (2011). In any case, the Lupus star-forming regions represent an excellent test-bed for 



theories of circumstellar disk evolution as their stellar members should span all evolutionary 
stages. 



The Scorpius clouds (Nozawa et al. 1991 Vilas-Boas et al. 2000) lie on the edge of 



the Lupus-Scorpius-Centaurus OB association, just north of the well studied Ophiuchus 
molecular cloud, but it is highly fragmentary and presents much lower levels of star- 
formation. In fact, the Gould Belt project only finds 10 YSOs candidates in the 2.1 sq deg. 
mapped by IRAC and MIPS (Hatchell et al., in preparation). The age of Sep is estimated 



to be ~ 5 Myr (Preibisch et al. 2002). 



The CrA star- forming region, also mapped by the Gould Belt project (Peterson et al 



2011), contains an embedded association known as the Coronet, a relatively isolated cluster 



containing HAeBe stars and TTauri stars (Chen et al. 1997). It is situated at a distance 



of 150 ± 20 pc out of the Galactic plane, at the edge of the Gould Belt (see Sicilia-Aguilar 



et al. 2008, and references therein). With an age of ~ 1 Myr, the Coronet is younger than 
the Lupus clouds and has been claimed to host an intriguingly high fraction of classical 



transition disks of ~ 50% (Sicilia-Aguilar et al. 2008). However, Ercolano et al. (2009) 



convincingly shows that the dust emission in T Tauri stars of spectral type M is very 



small short ward of 6/mi which might mimic an inner hole, the defining feature of typical 



transition disk systems. Based on this finding, Ercolano et al. (2009) estimate a much 



smaller fraction of transition disks in Coronet, of ~ 15%. 

2.1. Target selection 

We have systematically searched the catalogs of the c2d and Gould Belt Legacy 
Project^ applying the broad transition disk definition described in detail in Paper I to the 
Lupus I, III, IV, V, VI, Sep, and CrA clouds. In brief, we select systems that fulfill the 
following criteria. 

1. Have Spitzer colors [3.6]-[4.5] < 0.25, which excludes "fulls disks", i.e., optically thick 
disks extending inward to the dust sublimation radius except in cases with significant 



dust settling in inner disks around M stars (Ercolano et al. 2009). 



2. Have Spitzer colors [3.6]- [24] > 1.5, to ensure that all targets have very significant 
excesses (> 5-10 a), unambiguously indicating the presence of circumstellar material. 

3. Have S/N > 7 in 2MASS, IRAC, MIPS (24 //m) bands to only include targets with 
reliable photometry. 

4. Have K s < 11 mag, driven by the sensitivity of our near-IR Adaptive Optics 
observations and to avoid extragalactic contamination. 



5. Are brighter than R = 18 mag according to the USNO-B1 (Monet et al. 2003), driven 



by the sensitivity of our optical spectroscopy observations. Compared to the c2d 



2 the former is available at 



http://irsa.ipac.caltech.edu/data/SPITZER/C2D/ 



sample discussed in Merin et al. (2008) our sample might be slightly biased against 



very low mass stars and deeply embedded objects because of this brightness limit. 

These selection criteria result in a primary target list of 60 objects that we did follow-up 
using different observational facilities to characterize our transition disk candidate sample. 

3. Observations 

We performed multiwavelength (optical, infrared, and submillimeter) observations of 
our targets with the aim to identify which physical process is primarily responsible for their 
transition disk nature. High resolution optical spectra can be used to estimate spectral 
types and accretion rates from the velocity dispersion of the Ha. Near-IR images allow to 
identify multiple star systems down to projected separations of 0.06-007", corresponding 
to 8 — 14 AU at distances of 130 — 200 pc. From single dish submillimeter observations 
we inferred disk masses. In the following Section, we describe in detail the observations 
performed and the data reduction. 

3.1. Optical Spectroscopy 

We obtained high resolution (R > 20,000) spectra for our entire sample using 2 different 
telescopes: Magellan/Clay and Du Pont located at Las Campanas Observatory in Chile. 

3.1.1. Clay-Mike Observations 

We observed 49 of our 60 targets with the Magellan Inamori Kyocera Echelle (MIKE) 
spectrograph on the 6.5-m Clay telescope. The observations were performed on 2009 
April 27-28 and 2010 June 11-13. Since the CCD of MIKE's red arm has a pixel scale of 
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0.13"/pixel, we binned the detector by a factor of 3 in the dispersion direction and a factor 
of 2 in the spatial direction, thus reducing the readout time and readout noise. We used an 
1" slit width. The resulting spectra covered 4900-5000A at a resolving power of 22,000. 
This corresponds to a resolution of ~0.3 A at the location of the Ha line, and to a velocity 
dispersion of ~14 km s _1 . 

For each object, we obtained a set of 3 or 4 spectra, with exposure times ranging from 
3 to 10 minutes each, depending on the brightness of the targets. The data analysis was 
carried out with IRAFq] software. After bias subtraction and flat-field corrections with Milky 
Flats, the spectra were reduced using the standard IRAF package IMRED:ECHELLE. 

3.1.2. Du Pont-Echelle Observations 

The remaining 11 targets were observed with the Echelle Spectrograph on the 2.5-m 
Irenee du Pont telescope. The observations were performed in 2009 May 14-16, and we 
used an 1" slit width. The CCD's scale is 0.26" /pixel, and we consequently applied a 2x2 
binning. The wavelength coverage of the obtained spectra ranged between 4000 and 9000 A 
at a resolving power of 32,000 in the red arm. This corresponds to a resolution of ~0.2A 
and a velocity dispersion of ~9.4 km s _1 in the vicinity of Ha. 

For each object we obtained a set of 3 to 4 spectra with exposure times ranging from 
10 to 15 minutes each, depending on the brightness of the target. The data analysis was 
carried out with IRAF. After bias subtraction and flat-field corrections with Milky Flats, 
the spectra were reduced using the standard package IMRED:ECHELLE. 

3 Image Reduction and Analysis Facility, distributed by NOAO, operated by AURA, Inc., 
under agreement with NSF 
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3.2. Adaptive Optics Imaging 

High spatial resolution near-IR observations of our 60 targets were obtained with NaCo 
(the Nasmyth Adaptive Optics Systems (NAOS) and the Near-IR Imager and Spectrograph 
(CONICA) camera at the 8.2-m telescope Yepun), which is part of the European Southern 
Observatory's (ESO) Very Large Telescope (VLT) in Cerro Paranal, Chile. The data were 
acquired in service mode during the ESO's observing period 083 (2009 April 1 - September 
30). 

To take advantage of the near-IR brightness of our targets, we used the infrared 
wavefront sensor and the N90C10 dichroic to direct 90% of the near-IR light to the adaptive 
optics systems and 10% of the light to the science camera. We used the S13 camera 
(13.3 mas/pixel and 14x14" field of view) and the Double RdRstRd readout mode. The 
observations were performed through the K s and J-band filters at 5 dithered positions per 
filter. The total exposure times ranged from 1 to 50 s for the K s -band observations and 
from 2 to 200 s for the J-band observations, depending on the brightness of the target. The 
data were reduced using the Jitter software, which is part of ESO's data reduction package 
Eclipse^ 

3.3. Submillimeter Wavelength Photometry 

As discussed in the following section, our spectroscopic observations showed that our 
initial sample of 60 transition disk candidates was highly contaminated by asymptotic 
giant branch (AGB) stars. The 17 bona fide PMS stars were observed with the Atacama 



http : //www . eso . org/projects/aot/eclipse/ 
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Pathfinder Experiment (APEX) FJ the 12-m radio telescope located in Llano de Chajnantor 
in Chile. The observations were performed during period 083 (083.F-0162A-2009, 9.2 hrs) 
and period 085 (E-085.C-0571D-2010, 30.9 hrs). We used the APEX-LABOCA camera 



( Siringo et al.| 2009) at 870 /im (345 GHz) in service mode aiming for detections of the 
dust continuum emission. The nominal LABOCA beam is full width at half-maximum 
18.6 ± 1.0" and the pointing uncertainty is 4". To obtain the lowest possible flux limit, 
the most sensitive part of the array was centered on each source. The observations were 
reduced using the Bolometer array data Analysis package Boi 4 ^ 



For both observing runs, Skydips were performed hourly and combined with radiometer 
readings to obtain accurate opacity estimates. The absolute flux calibration follows the 



method outlined by Siringo et al. (2009) and is expected to be accurate to within 10%. 
The absolute flux scale pointing calibrators were determined through observations of either 
IRAS16342-38 or G34.3 while planets were used to focus the telescope. The telescope 
pointing was checked regularly with scans on nearby bright sources and was found to be 
stable within 3" (rms). 

The period 083 observations were performed using compact mapping mode with raster 
spiral patterns. The weather conditions were excellent with precipitable water vapor levels 
below 0.5 mm. Eight sources (objects # 1, 2, 5, 9, 12, 15, 16, 17) were observed. On-source 
integrations of 64 minutes were performed to achieve an rms ~ 7 mJy/beam. The brightest 
object of the whole sample (# 2) was the only source detected at submillimeter wavelengths 
in period 083. During the longer period 085 observing run, the beam switching mode 



5 This publication is based on data acquired with APEX which is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern Observatory, and the 
Onsala Space Observatory. 

6 http://www.apex-telescope.org/bolometer/laboca/boa/ 
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using the wobbling secondary and mapping mode were employed. During this period, 
the remaining nine sources were observed and object #12 was re-observed with higher 
sensitivity. The weather conditions were favorable with precipitable water vapour levels 
below 1.2 mm. The wobbler observations of each target consist of a set of two loops of 10 
scans per target, reaching a total on-source observing time of 48 minutes. An average rms 
~ 4 mJy/beam was obtained. In the case of a signal detection on-source position, we took 
a few maps in order to check for emission contamination from the off-position. In all cases, 
the contamination was discarded and we confirmed the detection of six sources (# 3, 7, 8, 
10, 11, and 12). 

4. Results 

4.1. AGB Contamination 

AGB stars are surrounded by shells of dust and thus have small, but detectable, 
IR excesses. The Spitzer-selected YSO samples from c2d and Gould Belt catalogs are 
therefore contaminated by AGB stars. Using high resolution optical spectra, we discovered 
that 43 objects of our candidates are AGB stars, while the remaining 17 targets are 
spectroscopically confirmed TTauri stars. We separated contaminating AGB stars from 
genuine transition disk TTauri stars in the same way as in Paper I, i.e. based on the 
presence/absence of emission lines associated with chromospheric activity and/or accretion 
and the presence of the Li 6707 A absorption line indicating stellar youth. The coordinates, 
Spitzer names, the USNO-B1 R-band magnitude, and the near to mid-IR fluxes of the AGB 
stars contaminating our sample of transition disks are compiled in Table [TJ 

As shown in Table [2] the fractional contamination due to AGB stars of our color 
selected transition disk candidates differs significantly between the different clouds. The 
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number of transition disk candidates is far too small in the case of Lupus I and Sep to draw 
any conclusions. Our Lupus III, IV, and CrA samples are contaminated by a fraction of 
AGB stars that is more or less consistent with the contamination in Ophiuchus (see Paper I, 
section 4.1.2). 

The small number of transition disks in CrA seems to be in contradiction with the 
larger sample identified by Sicilia-Aguilar et al. ( |2008 ). However, our selection criteria 
contain relatively strong brightness constraints (in particular K < 11) due to the design 
of our follow-up program which excludes most of the systems listed by them. In addition, 
as mentioned in the introduction, a large fraction of the transition disk candidates of 



Sicilia-Aguilar et al. (2008) might be classical M-dwarf TTauri stars with intrinsically little 



near IR excess due to the small color contrast between the disk and the stellar photosphere 



(Ercolanoet al. 2009) 



Apparently, Lupus V and VI are dramatically more contaminated than Lupus III, 
IV, i.e. all the color-selected transition disk candidates are in fact AGB stars. This high 
percentage of contamination is perhaps related to the position in the Galaxy (see Table [2]). 
The Lupus complex occupies 334 < 1 < 352, +5 < b < + 25, i.e. observing Lupus V and 
VI we are looking towards the galactic center closer to the plane. In contrast, CrA and 
Ophiuchus (Paper I) are located at higher Galactic latitudes. In any case, the absence of 
any spectroscopically confirmed transition disk in Lupus V, VI puts doubts on the finding 



of Spezzi et al. (2011, see their section 5.1) that the high fraction of Class III Lada systems 



can not be explained by contamination. So far all Class III YSO candidates from these 
clouds that have been followed up spectroscopically are clearly contaminating background 
giants. Our sample of transition disks in Lupus V, VI shares 30 Class III objects and one 



Class II object with the sample investigated by Spezzi et al. (2011). All these 31 objects 



turned out to be AGB stars which means that at least ~ 50% and potentially much more 
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of the Class III objects from Spezzi et al. (2011) are not YSOs but giant stars. This result 



also questions the conclusion of Spezzi et al. (2011) that Lupus V, VI are significantly older 
than Lupus I, III. 

4.2. Color selection of AGB star candidates 

The generally large fraction of giant stars in our sample of transition disk candidates 
allows to investigate possible refinements of our color selection algorithm. Figure [I] shows 
the color-color diagram of the 60 selected southern transition disk targets of this paper. 
All transition disk candidates in our sample with [3.6] — [24] < 1.8 turned out to be giant 
stars. This agrees quite well with the results obtained for the Ophiuchus sample where 
4/6 transition disk candidates with [3.6]- [24] < 1.8 had to be classified as giant stars (see 
Paper I, Fig. 1). Consequently, one may derive an estimate of the contamination of YSO 
catalogs due to background giant stars by applying this simple color cut. Figure[2] and [3] 
show the transition disk candidates and AGB candidates for both the c2d and GB catalogs 
including all star forming regions. Note, that we here apply color selection criteria only, 
i.e. the requirements of R < 18 mag and K < 11 mag that have been used to define 
the transition disk candidate sample for our multiwavelength follow-up program are not 
incorporated. Instead, here we are interested in estimating the fraction of YSO candidates 
in a given cloud that are likely to be AGB stars based on their very low 24/im excess. 

The resulting rough estimates of giant star contamination are given in Table [3] separated 
by catalog and cloud. According to these estimates, the AGB contamination is expected 
to vary significantly ranging from ~ 1 — 85 per cent. This shows that AGB contamination 
can have an important impact on studies that are based on the pure numbers of YSOs as 
provided by the c2d and GB catalogs. For example, star formation rates as determined 



e.g. in Heiderman et al. (2010) might become significantly smaller if AGB contamination is 
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taken into account. 

Apparently, applying the new more restrictive color selection could also significantly 
increase the success-rate of identifying YSOs directly from color selection criteria and future 
follow-up studies may take this into account. 

4.3. Characterizing southern transition disks 

The Spitzer and alternative names, 2MASS and Spitzer fluxes, and the USNO-B1 
R-band magnitudes and the relevant information derived from our follow-up observations 
for the remaining 17 bona fide transition disk candidates are listed in Tables [4] and [5} In 
what follows, we use the data discussed in Section [3] to characterize our sample of transition 
disks. 

4-3.1. Spectral types 



In oreder to determine the spectral types of the transition disks in our sample we use 



the equations by Cruz h Reid (2002) that empirically relate the spectral type with the 



strength of the Ti05 molecular band. The uncertainty of this method is estimated to be 
~0.5 subclasses. For most of our transition disk objects, estimates of the spectral types 



have been provided previously (Hughes et al. 1993, 1994; Krautter et al. 1997 Walter et al. 



1997 Sicilia-Aguilar et al.||2008| . The spectral types obtained by us and those given in the 
literature are listed in Tableland we find good agreement. All but one system (# 2) have 



been classified as M-dwarfs. For target # 2, we adopt the spectral type K0 given by Hughes 



eTaL1(fl993l). 



17 



4-3.2. Multiplicity 



Binarity can play an important role in the context of transition disks as the presence of 
a close stellar companion may cause the inner hole, i.e. some of the transition disks in our 
sample might actually be nothing else but circumbinary disks. Some systems in our sample 



have been previously identified as wide binaries. Merin et al. (2008) carried out an optical 
survey of the Lupus I, III, IV regions using bands Re, Ic, ZwI of the Wide-Field Imager 
(WFI) attached to the ESO 2.2-m telescope at La Silla Observatory. Visual inspection of 
the images revealed that objects # 1, # 3, # 4, and # 9 are wide binary systems with 
companions at projected separations of 420, 1000, 600, and 560 AU considering the distance 
of the Lupus clouds. 

We have newly identified six multiple systems by visual inspection of the NaCo images, 
objects # 6, # 7, # 11, # 13, # 16, and # 17 (see Figure [5]). The projected separations 
are 0.7", 0.4", 1.15", 1.8", 0.5", 0.5" corresponding to 140, 80, 230, 234, 75, 75 AU at the 
corresponding distances. Object # 13 is a triple system, i.e. a binary with an additional 
faint companion at 3" (390 AU). 

For each binary system, we searched for additional tight companions by comparing 
each other's point-spread functions (PSFs). The PSF pairs were virtually identical in all 
cases, except for target #11. The south-west component of this target has a perfectly 
round PSF, while the south-east component, ~1.5" away, is clearly elongated (see Figure 
111). Since variations in the PSF shape are not expected within such small angular distances 
and this behavior is seen in both the J- and K s -band images, we conclude that target #11 
is likely to be a triple system. 

We have also searched in the literature for additional companions in our sample that 
our VLT observations could have missed. In addition to multiple systems discussed above, 



we find that # 14 has been reported by Kohler et al. (2008) as a binary system with a 



projected separation of 0.13" (corresponding to 20 AU) and flux ratio of 0.7 using speckle 
interferometry at the New Technology Telescope in 2001. We see no evidence for a bright 
companion in our NaCo images (see Figure 111). However, since our AO images were taken 8 
years later than the speckle data, it is possible that the projected separation had changed 
enough in the intervening years for the binary to become unresolved. 

Hence, our sample consists of eleven multiple systems, i.e. nine binaries (objects # 1, 
3, 4, 6, 7, 9, 14, 16, 17), and two triples (objects # 11 and 13). Only in the cases of the close 
B/C pair in object #11 and # 14 the binary separation is small enough to suspect that 
the companions might have tidally disrupted the circumbinary disk thereby causing the 
inner hole inferred from the SED. However, in neither case the circumbinary nature can be 
confirmed because it is unknown whether the IR excess in object #11 originates in the wide 
A component or the close B/C pair and the multiplicity of object # 14 is not confirmed 
by our observations. We therefore only consider these two objects to be circumbinary disk 
candidates. Table [5] lists the projected angular separations of the systems. 

4.3.3. Disk Masses 



Andrews & Williams (2005, 2007) modeled the IR and submillimeter SEDs of 



circumstellar disks and found a linear relation between the submillimeter flux and the disk 



mass that has been calibrated by Cieza et al. (2008) who obtained 



o ,.^(0.86 mm) , rfi N9n 

M DISK ~ 8.0 x 10~ 2 [-^—^ L x — ?— 2 M JUP , 1 

mJy 140 pc 

where d± is the distance to the target. As described in Paper I, disk masses obtained with 
the above relation are within a factor of 2 of model derived values, which is certainly good 
enough for the purposes of our survey. However, one should keep in mind that larger 



systematic errors can not be ruled out (Andrews & Williams 2007) as long as strong 
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observational constraints on the grain size distributions and the gas-to-dust ratios are 
lacking. 



Adopting distances of 150 pc to Lupus IV, 200 pc to Lupus III ( Comeron||2008 ), 130 pc 



to Sep (Hatchell et al. in preparation), and 150 pc to CrA ( Sicilia-Aguilar et al.|2008 ) we use 



Equation ([!]) to estimate disk masses for the 17 systems in our sample (see Table [6]). 50% 
(i.e. 7/17) of the targets have been detected at 8510 jum (Table [5]). The corresponding disks 
masses range from 1 — 9 Mjyp. Adopting a flux value of 3<r for targets with non-detected 
emission, we derive upper limits for the remaining targets of ~ 1 — 4 Mjup. Most of our 
targets have disk masses < 1 — 2 Mjup, but 5 targets have disk masses typical for CTTSs 
(~3 - 10 Mjup). The most massive disks are detected around objects # 2 and # 3, with 9 
and 6 Mjup, respectively. 

4-3.4- Accretion rates 

The accretion rate is the second crucial parameter necessary to distinguish between 
the different mechanisms that may form inner opacity holes in circumstellar transitions 
disks. Most PMS stars show Ha emission, either generated from chromospheric activity 



or magnetospheric accretion (Natta et al. 2004). While non-accreting objects show 



rather narrow (< 200 km s _1 ) and symmetric line profiles of chromospheric origin, the 
large- velocity magnetospheric accretion columns produce broad (> 200 km s _1 ) and 
asymmetric line profiles. As in Paper I we estimate the accretion rates of our transition disk 



systems according to the empirical relation obtained by Natta et al. (2004), i.e. 



log(M acc (M Q jt- 1 )) = -12.89(±0.3) + 9.7(±0.7) x 10" 3 AV(km s" 1 ) (2) 

which is supposed to be relatively well calibrated for velocity widths covering 
AV = 200 — 600 km s -1 , which corresponds to accretion rates of 10~ n M yr _1 < 
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M acc < 10~ 7 M yr^ 1 . However, the empirical dividing line between accreting and 
non-accreting objects has been placed slightly shifted by different authors at A V between 



200 km s _1 ( Jayawardhana et al. 2003) and 270 km s _1 (White & Basri 2003). For systems 
with A V ~ 200 — 300 km s _1 we therefore separate accreting and non-accreting objects 
based on the (a) symmetry of the Ha emission line profile and take into account the spectral 
type because accreting lower mass stars tend to have narrower Ha emission lines. 

To measure the Ha velocity width AV we considered for each system the spectral 
range that corresponds to Ha ± 2500 Km/s. The continuum plus emission profile were 
fitted with a Gaussian plus parabolic profile. The parabolic fit was then used to normalized 
the spectrum. A single Gaussian profile was sufficient here, being the emission single or 
double-peaked, as at this stage we were only interested in obtaining a good parabolic fit for 
the normalization. Once the continuum had been normalized we measured A V at 10 per 
cent of the peak intensity. The obtained velocity dispersion of the Ha emission lines and 
the corresponding accretion rate estimates are given in Table [5] and Table [6j respectively. 
The obtained accretion rates should be considered order-of-magnitude estimates due to the 
large uncertainties associated with Equation (pi) and the intrinsic variability of accretion in 
T Tauri stars. 

Our sample shows a large diversity of Ha signatures. Five targets are classified as non- 
accreting objects that clearly show symmetric and narrow Ha emission line profiles (< 200 
km s _1 , see Figure p3J) as expected from chromospheric activity. For all these non-accreting 
objects, we estimate an upper limit of the accretion rate, i.e. M acc < 10~ n M yr _1 . 

We classify the remaining 12 transition disk objects as accreting. The majority of them 
(10) show clearly broad and asymmetric emission-line profiles (see Figure |6J). However, 
targets # 8 and # 12 represent borderline cases with a rather small velocity dispersion 
for accreting systems (A7~ 200 km s^ 1 ) and not clearly asymmetric line profiles. Such 
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borderline systems require a more detailed discussion. Both stars are of late spectral type 
(M4-M5) and very low-mass stars tend to have narrower Ha lines than higher mass objects 



because of their lower accretion rates (Natta et al. 2004) and their lower gravitational 



potentials (Muzerolle et al. 2003). Object # 12 additionally shows 8/im excess emission 
indicating the presence of an inner disk. Given all the available data, we classify targets 
# 8 and # 12 as accreting objects, but warn the reader that their accreting nature is less 
certain than that of the rest of the objects classified as CTTSs. As accretion in T Tauri 
stars can well be episodic, mutli-epoch spectroscopy would be useful to unambiguously 
identify the accreting nature of these two systems. 

The mass accretion rates estimated for the 12 disks classified as accreting systems 
range from 10 -11 M yr _1 to 10~ 7 ' 7 M yr _1 

5. Discussion 

5.1. The origin of the inner opacity hole 

With the collected information presented in the previous sections we have at hand the 
following information of the Spitzer-selected transition disks in our sample: 



Detailed SEDs that we quantify with the two-parameter scheme introduced by |Cieza" 



et al. (2007), which is based on the longest wavelength at which the observed flux is 



dominated by the stellar photosphere, A tumo fln and the slope of the IR excess, a, 



To calculate A tU rnoff, we compare the extinction-corrected SED with NextGen Mod- 



els (Hauschildt et al. 1999) normalized to the J-band and choose A turno fj as the longest 
wavelength at which the stellar photosphere contributes over 50% of the total flux. The 
uncertainty of A turno ff is roughly one SED point. 
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computed as dlog(F\)/d\og(\) between A tU rnoff and 24 /xm. 

• Multiplicity information from the literature and AO IR imaging. 

• Disk mass estimates based on measured submillimeter flux. 

• Accretion rate estimates derived from Ha line profiles. 

This information allows to separate the sample according to the physical processes 
that are the most likely cause of the inner opacity hole: grain growth, planet formation, 
photoevaporation, or close binary interactions. In what follows we briefly review each 
process that might be responsible for the formation of the inner opacity holes, describe 
our criteria for classifying transition disks, and discuss the corresponding sub-samples of 
transition disks obtained. 

5.2. Accreting transition disks 

The presence of accretion in classical transition disk objects raises the question how 
the inner disk can be cleared of small grains while gas remains in the dust hole. The two 
mechanisms that can explain the coexistence of accretion and inner opacity holes are grain 
growth and dynamical interactions with (sub)stellar companions. 

5.2.1. Grain growth 

Due to both, the higher densities and the faster relative velocities of particles in the 
inner parts of the disk, this disk region offers much better conditions for dust agglomeration 
than the outer parts of the disk. Therefore, significant grain growth should start in the 
inner disk regions. As soon as the grains grow to sizes exceeding the considered wavelength 



23 



(r >> A), the opacity decreases until an inner opacity hole forms. Early models by 



Dullemond & Dominik (2005) taking into account only dust coagulation predict much too 
short timescales of the order of 10 5 yrs to clear the entire disk of small grains, which is 
inconsistent with observed SEDs of most classical T Tauri stars. A more reliable picture 
combines coagulation and collisional fragmentation or erosion of large dust aggregates 



(Dominik & Dullemond 2008). 



As a gradual transition between the inner and outer disk is predicted by grain growth 



and dust settling models (7Weidenschilling 2008), grain growth dominated disks should 



have a exceS s < (i.e., falling mid-IR SEDs) while A tur noff, associated with the size of the 
hole, can differ over a rather wide range of values. Although grain growth does not directly 
affect the gas, it may increase accretion because the inner opacity hole can lead to efficient 



ionization and trigger the MRI instability (Chiang & Murray-Clay 2007) 



Among the 17 transition disks in our sample, nine objects are accreting and are 
associated with a excess < 0. The corresponding SEDs are shown in Figure [7| The grain 
growth candidate systems in our sample could be confounded with accreting classical 



T Tauri M-stars as predicted by Ercolano et al. (2009). However, most of our grain growth 



dominated disks have SEDs close to the stellar photosphere up to ~ 8 — 10 /mi and 
we therefore do not expect significant contamination by non transition disks. The only 
exceptions being objects # 5 and 10 with a small value of Aturnoff— 2.2 and to some extend 
# 4 and #11 (A tur noff= 4.5). We recommend the reader to keep in mind the uncertain 
classification of these four systems. 

Compared with the Oph sample (Paper I) the accretion rates obtained for grain growth 
dominated disks are slightly smaller (i. e. 10~ 93 — 1O _11 M yr _1 ). This might be related to 



the slightly lower stellar masses or to advanced viscous evolution as discussed in |D'Alessio 
et~all([2006|. 
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The grain growth dominated disks are located in the Lupus III, IV (8) and the CrA 
4) star- forming regions. 

5.2.2. Dynamical interactions with (sub)stellar companions 
The truncation of the disk as the result of dynamical interactions with companions 



was first proposed by Lin & Papaloizou (1979). More recently, it has been shown that 



most PMS stars are in multiple systems with a lognormal semi-major axis distribution 



centered at ~ 30 AU (e.g., Ratzka et al. 2005). A significant fraction of the binaries in the 
star-forming regions considered here should therefore be tight binaries with separations of 
1 — 20 AU. Disks in such close binary system will be tidally truncated at ~ 2 x the binary 



separation and a circumbinary disk with an inner hole is formed (Artymowicz & Lubow 



1994). The corresponding SED is that of a transition disk. However, the circumbinary 



nature does not exclude additional evolutionary processes to be at work and we therefore 
provide an additional classification based on the disk properties only (see Table [6]). 



Identification of companions that may cause the formation of a circumbinary disk is 
possible either due to high-resolution imaging or by measuring radial velocity variations. As 



described in Sect. |3.2[ we identified 2 circumbinary disk candidates among our 17 transition 
disk systems. One of them, object # 11 was discovered by inspecting the NaCo images 
obtained with the VLT, while the close binary nature of object # 14 has been discovered 



by Kohler et al. (2008) using speckle interferometry at the NTT. 



Object #11 shows signs of strong accretion and has a SED with a excess ~ — 1 in 
agreement with grain growth. It is currently not clear under which conditions gap-crossing 
streams can exist and allow accretion onto the central star to proceed, but signs of accretion 



in circumbinary systems (Carr et al. 2001; Espaillat et al. 2007) indicate that accretion is 
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likely to continue. On the other hand, object # 14 is a non-accreting system such as the 



known binary CoKu Tau/4 (Ireland & Kraus 2008). According to its Ldisk/L* ratio we 



classify this system as a circumbinary/photoevaporation disk candidate. 

As a final note of caution, we would like to stress that both objects discussed above (# 
11 and # 14) are circumbinary disk candidates. As all but one of our targets are M-type 
stars, most companions potentially responsible for their transition disk SEDs are expected 
to lie at closer separations than those probed by the AO images. Therefore, our sample of 
circumbinary disk candidates is incomplete and heavily biased towards large separations. 
Methods more sensitive to closer companions such as aperture masking and/or radial 
velocity observations are required to draw firm conclusions on circumbinary disk fractions. 

5.2.3. Giant Planet formation 

The most exciting way to produce a transition disk SED is by giant planet formation. 
According to early models as well as recent numerical simulations, the formation of giant 



planets involves the formation of gaps and holes in the circumstellar disk (Lin & Papaloizou 



1979 Artymowicz & Lubow 1994) 



As in the case of (sub) stellar companions it is uncertain if and to what extent 
accretion proceeds in the presence of a forming giant planet. Therefore, the most important 
sign of ongoing planet formation remains a sharp inner hole, usually corresponding to 
«excess> (i.e., a rising mid-IR SED). However, although very useful, the definition of 
^excess is incomplete, as the SED may also steeply rise at wavelengths longer than 24 //m. 
A spectacular example illustrating this is given by object # 3. While a ex cess~ — 2.2, 
the SED steeply rises between 24 /j,m and 70 /im. Furthermore, object #3 shows clear 
signs of accretion (M acc = 10~ 10,1 M /yr) and of harboring a relatively massive disk 
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(M DISK ~ 6 Mjup). Since this is a very atypical object, we verified that the large 70 /im flux 
is not contaminated by extended emission from the molecular cloud. We examined the 24 
and 70/zm mosaics and verified that the detections are consistent with point sources at the 
source location (see Figure |8|. A more typical transition disk system that might represent 
a currently planet forming disk is target # 15 with a clearly positive value of a ex cess an d 
a high accretion rate. A borderline case between grain growth and planet forming disks 
is object # 2. A high accretion rate combined with a exC ess~ could be consistent with 
both scenarios. Keeping in mind the ambiguity, we classify object # 2 asa planet forming 
disk candidate because it could potentially be an extremely interesting object. The SED 
of object # 2 might be explained by a discontinuity in the grain size distribution rather 
than an inner opacity hole. While the inner part of the disk still contains small grains, 
outer regions of the disk might be dominated by slightly larger dust particles. Such a 
scenario is in excellent agreement with the predictions of numerical simulations performed 



by Rice et al. (2006). They show that the planet-disk interaction at the outer edge of the 
gap cleared by a planet can act as a filter passed by small particles only which produces a 
discontinuity in the dust particle size. To firmly establish its nature object # 2 deserves 
further follow-up observations (e.g., high resolution imaging with ALMA). 

The SEDs of the three candidates for ongoing giant planet formation in our sample are 
shown in Figure [9] The hosting forming giant planets candidates are located in the Lupus 
III, IV (2) and the CrA (1) star-forming regions. 



5.3. Non-accreting objects 

The second main class of transition disks are those that do not show signs of accretion. 
In such disks the inner opacity hole, i.e. the lack of small dust particles in the inner disk 
regions, is likely to coincide with a gas hole, i.e. the inner disk is completely drained. 
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5.3.1. Photoevaporating disk 



The most important process for clearing the inner disk in transitions disks that do 



not accrete is photoevaporation (e.g. Alexander et al. 2006). According to this model 



extreme-ultraviolet (EUV) photons, originating in the stellar chromosphere, ionize and heat 
the circumstellar hydrogen which is then partly lost in a wind. This process is supposed 
to work in all circumstellar disks but becomes important only when the accretion rate 
drops to values similar to the photoevaporation rate. Then, the inner disk drains on the 



viscous timescale supported by the generation of the MRI (Chiang & Murray-Clay 2007). 



Once an inner hole has formed, the inner disk rim is efficiently radiated and the entire disk 
should therefore quickly disappear. Photoevaporating disks should have negligible accretion 



(Williams & Cieza 2011). To separate photoevaporating disks from debris disks, we require 



the disk luminosity to be higher than that seen in the brightest bona-fide debris disks, 



i.e. L DISK /L* > 10~ 3 ( |Bryden et aT1|2006| |Wyatt||2008[ ). We thus obtained 70 fjcm upper 
limits from the noise of the 70 /zm Spitzer mosaics at the source location and calculated 
Ldisk/L* for our sample by integrating the stellar fluxes and disk fluxes over frequency (see 
Section 5.1.3 in Paper I for details of the 70 /zm data analysis and the L D isk/L* calculation). 

We classify three transition disks as photoevaporating disk candidates with 
negligible accretion (M acc < 10~ n yr _1 ) but Ldisk/L* > 10~ 3 (Table ^j). According to our 
submillimeter measurements, all these three systems have small disk masses (< 1 — 3.4 Mjup, 
Table [6]). In fact, for all photoevaporation candidates we could only derive upper limits on 
the disk mass. The SEDs of the three systems classified as photoevaporating disk objects 



are shown in Figure 10 The photoevaporated disks are located in the Lupus III, IV (2) and 



CrA (1) star- forming regions. 



-28- 

5.3.2. Debris disk 

Photoevaporation can be considered as a transition stage between primordial and 
debris disks. Debris disks contain a small amount of dust and are gas-poor. We find two 
debris disk candidates, i.e. non-accreting systems with Ldisk/L*< 10~ 3 , among our 17 



transition disks (see Figure 11). The debris disks are located in the CrA (1) and Sep (1) 
star-forming regions. 

5.4. Implications for disk evolution 

5.4-1. Heterogeneity of transition disks 

In the previous sections, we presented detailed follow-up observations of 60 Spitzer- 
selected transition disk candidates located in the southern star-forming regions Lupus I, 
III, IV, V, VI, CrA, and Sep. Optical spectroscopy revealed that only 17 systems of these 
candidates are genuine transition disk T Tauri stars. Deriving estimates for the accretion 
rates, disk masses, and multiplicity of these 17 systems we classified them as dominated by 
grain growth (9), giant planet formation (3), photoevaporation (3), or being in the final 
debris disk stage (2). Two of these transition objects, one grain growth ($ 11) and one 
photoevaporating (# 14), are circumbinary disk candidates, which offers the possibility of 
tidal truncation as mechanism responsible for an inner hole in the common/shared disk. 
Combining these results with those presented in Paper I, we now have at hand well-defined 
and well characterized samples of transition disks from several different star-forming regions. 



Figure [12] summarizes the properties of these samples based on a exC ess and Atumoff- The 
main aim of these series of papers is to progress with our understanding of circumstellar 
disk evolution and to compare transition disk samples of different clouds is key in this 
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respect. Table [7] shows the fractions of different types of transition diskajfor Ophiuchus 



age ~ 0.3 — 2.1 Myr, Wilking et al. 2005, and references therein), CrA (age ~ 1 Myr, 



Sicilia-Aguilar et al. 2008), and Lupus I, III, IV (age 1.5 - 4 Myr, Merin et al. 2008). All 



YSO candidates followed up spectroscopically located in Lupus V, VI turned out to be 



AGB stars. These clouds have been recently estimated to be a few Myrs older (Spezzi et al. 



2011) based on the dominance of Class III systems. As we have shown in Section 4.1, at 
least ~ 50 % of the claimed class III systems located in Lupus V, VI are very likely to be 
AGB stars. This reduces the fraction of class III objects to values similar to those obtained 
for Lupus III. Based on this, Lupus V, VI, and III could well be of a very similar age. 

The main result that can be obtained from Table [7] clearly is that young clouds 
(< 1 — 4 Myr) contain a mixture of grain growth, photoevaporating, debris, and tidally 
disrupted transition disks. It is clear that all states of disk evolution are already present at 
this age range, which implies that different disks evolve at different rates and/or through 
different evolutionary paths. An important difficulty in constraining disk evolution is that 
stellar ages obtained from isochrones are very uncertain for individual systems. An analysis 
of the stellar age distributions of each disk category is therefore postponed to paper III 
(Cieza et al. 2012, ApJ submitted), where we discuss a larger sample of well characterized 
transition disk objects including the systems presented here. 

5.4-2. Evidence for low photoevaporation rates 

The general picture of photoevaporation is the following. In very young disks, the 
accretion rate largely exceeds the evaporation rate and the disk evolves virtually unaffected 



8 Circumbinary disks are included twice in the table as binarity does not exclude a second 
process to cause the inner opacity hole in the disk. 
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by photoevaporation. As the accretion rate is decreasing with time, the disk necessarily 
reaches the time when the accretion rate equals the photoevaporation rate and the outer 
disk is no longer able to resupply the inner disk with material. At this point, the inner disk 
drains on the viscous timescale (< 10 5 yr) and an inner hole of a few AU in radius is formed 
in the disk. The inner disk edge is now directly exposed to the EUV radiation and the disk 
rapidly photoevaporates from the inside out. 

Early models of EUV photoevaporation predict evaporation rates of 1CT 10 — 10 -9 M /yr 



(Hollenbach et al. 1994). More recent simulations taking into account X-ray (Owen et al. 



2010) and/or far-ultraviolet (FUV) irradiation (Gorti & Hollenbach 2009) in addition to 
the EUV photons, largely exceed these early predictions, reaching photoevaporation rates 



of the order of 10 8 M /yr (see also Armitage 2011). 



As in steady state accretion disks the mass transfer through the disk is roughly 
proportional to the mass accretion rate onto the star, a crucial prediction of the 
photoevaporation model is that high photoevaporation rates imply high disk masses at the 
time the inner disk is drained. In particular, models with efficient X-ray photoevaporation 
predict a significant population of relatively massive (~7 Mjup) non-accreting transition 



disks (Owenet al. 2011). 



Figure 13 shows the upper limits (derived from submillimeter non detections) on 



the disk masses of the photoevaporating transition disks in all the clouds we considered 
so far. Even taking into account uncertainties in our classification of photoevaporation 
candidates, it is evident that large numbers of non-accreting but massive disks do not 
exist. This indicates that photoevaporation is less efficient than predicted by the models 
described above. However, one has to take into account that the sample of transition disks 
considered here contains low-mass stars only while model calculations have been performed 
exclusively for more massive stars ~ 1M Q . Therefore, either a more homogeneous sample of 
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photoevaporating disk systems covering a larger range of host star masses (earlier spectral 
types) or simulations of photoevaporation for disks around low-mass stars are required to 
provide a final answer on this issue. 

5.5. Current limitations and future perspectives 

Of course, our classification of transition disk objects is based on rather rough empirical 
relations and requires to carefully consider possible caveats. An obvious uncertainty 
concerns our multiplicity survey. The method of direct detection of companions is 
obviously more sensitive to binaries with large separations and low inclinations. Our NaCo 
observations are sensitive to projected separations of ~ 10 — 15 AU given the distance to 
our targets and - depending on the intrinsic distribution of orbital separations - we may 
therefore miss a significant fraction of close binaries. To overcome this observational bias 
we are currently performing radial velocity measurements of our targets using VLT/UVES. 
The method of detecting radial velocity variations is more sensitive to small separations 
and high inclinations and therefore complements the imaging results presented here. We 
will present the results in a forthcoming paper. However, the fact that only 6 of the 43 
transition disks studied herein and in Paper I are circumbinary disk candidates strongly 
suggests that binaries at the peak of their separation distribution (~ 30 AU) do not result 
in transition disk objects as such stellar binaries would be easily detectable by our AO 



observations. Instead, they are likely to destroy the disk rather quickly (Cieza et al. 2009). 



Another uncertainty in our classification procedure is the rather ad-hoc separation 
between photoevaporating and debris disk systems by using a limit in Ldisk/L*. However, 
there is a physical and not only phenomenological difference between these two types of 
transition disks. Photoevaporating disks are dissipating primordial disks and should have 
gas rich outer disks while the debris disks should be gas poor. Molecular line observations 
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with the Atacama Large Millimeter/Submillimeter Array (ALMA) of non-accreting disks 
will be able to distinguish between the two types of objects. 

A huge problem related to the process of photoevaporation is that the mass loss 
rates predicted by different models differ by up-to two orders of magnitude (see e.g. 



Williams & Cieza 2011). The disk mass at the time photoevaporation opens a hole in 
the disk is directly connected to the photoevaporation rates. Measuring the disk masses 
of photoevaporating disks could therefore significantly constrain theoretical models of 
photoevaporation. However, all of the photoevaporting disk candidates remain undetected 
and we can only put upper limits to their masses. Fortunately, ALMA will be much more 
sensitive than all presently available telescopes and will soon be able to measure the masses 
of many bona fide photoevaporating disks. ALMA should also be able to measure, through 
high resolution continuum observations at multiple wavelengths, the radial dependence in 
the grain size distribution expected in the grain-growth dominated disks. 

Finally, the recent identification using the aperture masking technique of what seems to 



be forming planets within the inner cavities of the transition disks around T Cha (Huelamo 



et al. 2011) and LkCa 15 (Kraus & Ireland 2012) strongly encourages to obtain similar 



observations for the three planet-forming disk candidates identified herein, objects # 2, # 3, 
and # 15. Any system with a planet still embedded in a primordial disk would represent an 
invaluable laboratory to study planet formation with current and future instrumentation. 

6. Summary 

We have carried out a multifrequency study of Spitzer-selected YSO transition 
disk candidates located in the Lupus complex (53), CrA (5), and Sep (2). We obtained 
submillimeter observations (APEX), optical high resolution echelle spectroscopy (Clay/Mike, 
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Du Pont/echelle), and NIR images (from AO imaging VLT/NaCo). After deriving spectral 
types of each target, 43 AGB stars were removed (Lupus complex (41), CrA (1), and Sep 
(1)), leaving a sample of 17 genuine transition disk systems. We find that the vast majority 
of AGB stars have [3. 6] -[24] < 1.8, underscoring the need for a spectroscopic confirmation of 
YSO candidates with small 24//m excesses. The data obtained for the 17 transition systems 
allows to estimate multiplicity, stellar accretion rates, and disk masses thereby allowing 
to identify the physical mechanism that is most likely to be responsible for the formation 
of the inner opacity hole. The observational results of this study can be summarized as 
follows: 

1. The derived spectral classification indicates that all but one (object # 2, KO) central 



star are M-type stars, in agreement with previous results (Comeron 2008) 



2. 12/17 targets are accreting objects (i.e. asymmetric Ha profile having a velocity 
width > 200 km s _1 at 10% of peak intensity). 

3. ~ 50% of the sample are multiple systems and among them, two triple systems. Two 
binary systems have small projected separations and are therefore candidates to host 
a circumbinary disk. 

4. 7/17 targets have flux detection in the submillimeter. For the remaining systems, we 
derive and upper limit of the disk mass (corresponding to a flux of 3 x rms). The 
estimated disk masses for the detected objects cover the range 2 Mjup-10 Mjup. 

Combining the derived estimates of disk masses, accretion rates and multiplicity with 
the SED morphology and fractional disk luminosity (Ldisk/L*) allows to classify the disks 
as strong candidates for the following categories: 

• 9/17 grain growth-dominated disks (accreting objects with negative SED slopes in 
the mid-IR, a excess < 0). 



-34- 

• 3/17 photoevaporating disks (non-accreting objects with disk mass < 3 M JUP , but 
Ldisk/L*> 10- 3 ). 

• 2/17 debris disks (non-accreting objects with disk mass < 2.1 Mjup and 
Ldisk/L*< 10- 3 ). 

• 2/17 circumbinary disks (a binary tight enough to accommodate both components 
within the inner hole). 

• 3/17 giant planet forming disk (accreting systems with SEDs indicating sharp inner 
holes). 

Inspecting in more detail the different sub-clouds analyzed in this study we find the 
same heterogeneity of the transition disk population in Lupus III, IV, CrA as in our 



previous analysis of transition disks in Ophiuchus (Cieza et al. 2010, Paper I). We therefore 
conclude that photoevaporation, giant planet formation, and grain growth produce inner 
holes on similar timescales. Not a single transition disks has been found in Lupus I, V, VI. 
All 33 candidates that have been spectroscopically followed up turned out to be AGB stars 



which questions the recent interpretation of Spezzi et al. (2011) that Lupus I, V, VI might 



be relatively old star forming regions dominated by Class III objects. 

In addition, our detailed observational analysis of transition disks provides clear 
constraints on theoretical models of disk photoevaporation by the central star. According 



to the large evaporation rates predicted by recent models (i.e. see Armitage 2011), large 
numbers of massive photoevaporating transition disks systems should exist. In contrast to 
this prediction, all photoevaporating disk candidates identified in this work and Paper I 
contain very little mass, indicating much smaller evaporation rates at least for the low-mass 
stars considered here. Similarly, the low incidence of circumbinary transition disk candidates 
(~ 10%) supports the idea that most disks are destroyed rather quickly by companions at 
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~ 10 — 40 AU separations. 

Finally, we emphasize that the 43 transition disk systems discussed in this work and in 
Paper I represent the currently largest and most homogeneous sample of well-characterized 
transition disks. Further investigating these systems with new observing capabilities such 
as ALMA therefore holds the potential to significantly improvement our understanding of 
the physical processes driving circumstellar disk evolution. 
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Fig. 1.— Spitzer color-color diagram for the southern sample. Contaminating AGB stars 
greatly dominate the sample with [3.6] [24] < 1.8. 
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Fig. 2.— Spitzer c2d systems classified into YSO candidates, AGB candidates, and transi- 
tion disk candidates according to simple color cuts based on the results of our spectroscopic 
follow-up program (see text for more details). 
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Fig. 3— Same as Fig.[2]but for the GB catalog. 
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Fig. 4. — The K s -band images of the six multiple systems that have been detected with 
our VLT-AO observations (objects # 6, 7, 11, 13, 16, and 17) and of object # 14, which 
has been identified as a close binary with a 0.13" separation from speckle observations in 



2001 (Kohler et al. 2008). The putative companion in object # 14 remains unresolved by 



our 2009 observations. Targets #11 and 13 are triple systems. In the former case, the 
tighter components are not fully resolved, but their presence can be inferred from the highly 
elongated image (lower right panel). 
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Fig. 5. — Normalized average Ha profiles for non-accreting objects with clear Ha emission. 
The horizontal dashed line indicates the 10% peak intensity, where AV is measured. The 
velocity widths are < 200 km s _1 and the line profiles are symmetric. 
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Fig. 6. — Normalized Ha profiles of the 12 accreting objects. These systems are considered 
accreting because either the velocity width is > 200 km s" 1 or the line profile is asymmetric. 
Note that objects # 8 and 12 represent borderline systems as the AV~ 200km/s and the 
line is not clearly asymmetric (see text for more details). 
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Fig. 7. — SEDs of the 9 grain growth candidates. Disk masses range from < 1.1 to 5.1 Mjup, 
and accretion rates from 10~ n to 10 -9 ' 3 M yr _1 . The filled circles are detections, while the 
arrows represent 3<r limits. The open squares correspond to the observed optical and near-IR 
fluxes before being corrected for extinction as in Paper I. For each object, the average of 
the two R-band magnitudes (from the USNO-B1 catalog) listed in Tabled has been used. 
The classification of object # 5 and 10 as grain growth dominated is slightly uncertain as 



classical TTauris stars of late M-types can have similar SEDs (Ercolano et al. 2009). The 



solid line represents the stellar photosphere normalized to the extinction-corrected J band. 



The dotted lines correspond to the median mid-IR SED of CTTSs calculated by Fur Ian et al. 



(2006). The dashed lines are the quartiles. 
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Fig. 8. — The 24 and 70 /zm images of object # 3. We find no evidence for extended emission 
from the molecular cloud. In both mosaics, the detections are consistent with a point source 
at the location of the target (marked by the crosshairs). See the electronic edition of the 
Journal for a color version of this figure. 
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Fig. 9. — SEDs of the three planet forming candidates. While object # 15 can be considered 
a classical planet forming candidate system with awess > and A turno ff = 8.0 /mi, the other 
two systems are somewhat peculiar: target # 3 shows a very steep rise in flux observed 
at 70 /mi, which indicates a very large inner hole and object # 2 being relatively close to 
a full disk but with signs for a small and sharp inner hole. Disk masses are 9.1, 5.6 and 
< 2 Mjup for objects # 2, # 3, and # 15; respectively. The solid line as well as the dashes 
and dotted lines are the same as in Fig.[7j 
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Fig. 10. — SEDs of the 3 transition disks classified as photoevaporating disk candidates. The 
symbols, solid line as well as the dashes and dotted lines are the same as in Figure [7j None of 
these systems have been detected at submillimeter wavelength and consequently only upper 
limits for the disk masses could be derived. We conclude that photoevaporation seems to be 
less efficient than has recently been suggested. 
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Fig. 11. — SEDs of the 2 debris disk candidates identified in our sample. We distinguish 
between debris disks and photoevaporation applying an ad-hoc limit on Ldisk/L* ~ 10 -3 . 
The detection of the gas component in photoevaporating disks (e.g. with ALMA) may lead 
to a more physically motivated separation between the two sub-samples. The solid line as 
well as the dashes and dotted lines are the same as in Figure [7} 
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Fig. 12.— ^excess vs A tU rnoflf f° r transition disks identified in Paper I and the present work. The 
locations of the Ophiuchus transition disk sample from this work (top panel) and Paper I (bottom 
panel) cover very similar ranges in the a eX cess versus Atumoff plane. Different symbols indicate 
different formation processes of the inner opacity hole. In general, the Lupus, Cra, and Sep data 
confirm our previous findings: planet forming disks have large values of a ex cess and At urno ff, grain 
growth dominated disks should have small a exC css , but cover the entire range of Aturnoff ; and debris 
disks have extremely low values of a exC ess and the IR excess starts at long A turn off- However, two 
systems clearly show that a ex cess and A turno ff alone cannot fully characterize transition disks. We 
identified one planet forming disk candidate with a exceS s< but indications of a sharp hole at 
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Table 2. T Tauri and AGB star fractions 



Cloud 


Coordinates 


Sample 


Transition Disk 


AGB 


Percentage 


Age 


Cloud's 




(l,b) 




candidates 




of AGB 


[Myr] 


Distance 




~ (deg,deg) 


# 


# 


# 


% 




[pc] 


Lupus I 


339, 16 


2 





2 


100 


1.5-4 


150 ± 20 


Lupus III 


340, 9 


15 


10 


5 


33 


1.5-4 


200 ± 20 


Lupus IV 


336, 8 


5 


2 


3 


60 


1.5-4 


150 ± 20 


Lupus V 


342, 9 


16 





16 


100 


10 


150 ± 20 


Lupus VI 


342, 6 


15 





15 


100 


10 


150 ± 20 


Cra 


0,-19 


5 


1 


1 


20 


1 


150 ± 20 


Sep 


250,18 


2 


1 


1 


50 


5 


130 ± 20 


Oph (Paper I) 


353,18 


3d 


26 


8 


24 


2 


150 ± 20 
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Table 3: Distribution of YSOc, transition disks (TD) and AGBs candidates organized 
by c2d and Gould Belt Legacy Projects' clouds 



Region 


YSOc 


AGB 


TD 


AGB 






candidates 


candidates 


candidates in TD region 






whole sample 




in TD region 




# 


% 


# 


% 



c2d Legacy Project's clouds 



CHA II 


29 


LUP I 


20 


LUP III 


79 


LUP IV 


12 


OPH 


297 


PER 


387 


SER 


262 



10.4 
15 
18.9 
25 
7.7 
2.6 
6.5 



7 

8 
18 
5 

52 
56 
60 



Gould Belt Legacy Project's clouds 



28.6 

25 
11.1 

20 
15.38 
10.7 

10 



AURIGA 


174 


1.7 


28 


7.1 


CrA 


45 


4.4 


7 


14.2 


IC5146 


163 


2.4 


24 


4.1 


LUP V 


44 


47.7 


22 


36.3 


LUP VI 


46 


67.3 


21 


57.1 


SERP-AQUILA 


1442 


28.6 


641 


32.6 


CHAM I 


93 


1 


17 


5.9 


CHAM III 


4 


75 


1 


100 


MUSCA 


13 


84.6 


5 


80.0 


CEPH 


119 


2.5 


19 


10.5 


SCO 


9 


11 


4 


25 



0> 



go 

en 



a 
o 
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Table 5. Observed Properties 



# Ra (J2000) Dec (J2000) Tel. SpT 
(deg) (deg) 



SpT. References 3 - Li I b Ha b ' c Flux mm d crFlux mm Separ e 

(A) (km s -1 ) (mjy) (mJy) (arcsec) 



1 


240.10887 


-41.89877 


Clay 


M5.25,M1J 


1.2 


0.47 


162 


< 21 


7 


2.8 


2 


240.18554 


-41.92534 


Du Pont 


K0 


7 


0.44 


532 


100 


5 




3 


241.79833 


-39.06326 


DuPont 


M1.5,M0.5 


1.6 


0.41 


283 


34.5 


2.9 


5 


4 


241.96800 


-38.96840 


Du Pont 


M3.25.M1.5 


1.6 


0.46 


321 


< 9.9 


3.3 


3 


5 


242.05258 


-39.14264 


DuPont 


M2,M1.5 


1.6 


0.5 


233 


< 21 


7 




6 


242.11837 


-39.09229 


Du Pont 


M5,M4 


1.6 


0.28 


3-13 


< 10.8 


3.6 


0.7 


7 


242.13146 


-38.79148 


Du Pont 


M2.25,M2 


1,4 


0.25 


382 


6.7 


2.2 


0.1 


8 


242.17413 


-39.02695 


Du Pont 


M5.75,M5.5 


1.6 




200 


9.7 


2.5 




9 


242.23133 


-39.04276 


Du Pont 


M6,M6 


1.6 




189 


< 21 


7 


2.8 


10 


242.25583 


-39.41997 


Clay 


M4,M4 


1.2 


0.45 


369 


31.4 


3.4 




11 


242.47496 


-39.39109 


Du Pont 


M2.75,M1.5 


1.1 


0.40 


336 


16.7 


3.3 


1.15 


12 


242.62321 


-39.37076 


Clay 


M4.5,M4 


1.2 


0.52 


180 


23.2 


4.7 




13 


245.54000 


-19.88357 


Clay 


M3.7 


1 


0.55 


132 


< 10.8 


3.6 


1.8,3 


14 


285.12113 


-36.93437 


DuPont 


M4,M3 


1.5 


0.30 


93 


< 10.5 


3.5 


0.132 


15 


285.24187 


-36.75139 


Clay 


M0.75 


1 


0.48 


440 


< 21 


7 




16 


285.37088 


-37.03011 


DuPont 


M3.75,M3.5 


1.3 




83 


< 21 


7 


0.5 


17 


285.79929 


-37.15055 


Clay 


M5.0 


1 


0.51 


205 


< 21 


7 


0.5 



References. (1) Spectral type from this work; (2) |Merm et al.| ^2008) ; (3) |Sicilia-Aguilar et aT] | |2008| ; (4) |Krautter| 
|et al.| | |1997| l; (5) |WaIter et al.| l |1997| ;(6) |Hughes et al.| ( |1994| l; (7) |Hughes et al.| | |1993| l 

b " ■ ■ ■ " implies that the signal to noise in this region of the spectrum is too low to measure the width or establish 
the presence of the line 

c "-l" implies that Ha is seen in absorption. 

d "<" implies upper limits value (3 a) 

°Sources # 1, 3, 4, and 9 have been identified as binaries by WFI observations | |Merin et al. 2008). The binary nature of 
source # 14 was discovered using speckle intcrfcrometry at the NTT (La Silla) revealing a projected separation of 0.132" 
zt 0.009" jKohler et al.||2008| . Sources # 11, 13 are triple systems with tight binary components consistent with two 
equally bright objects and a projected separation of ~0.05". 

The spectral type derived by us is significantly later than the previously obtained value. We consider our estimate of 
M5.25 based on a high quality optical spectrum to be more reliable than the rather rough guess of|Merin et al.|(2008) 
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derived from optical and ncar-IR photometry only. 
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Table 6. Derived Properties 



# 


LOG(Acc. rate) 


Mass Disk 


Yproj. 


^■tun — off 


(^excess 


LOG(L D /L») 


L„ 


Classification 




(Mg 


/yr) 


(Mjc/p) 


(AU) 


/illl 






L Q 


Cloud partnership 


1 




< -11 


< 1.9 


120 


4.50 


-0.80 


-2.1 


0.06 


photo, disk, Lup IV 


2 




-7.7 


9.1 




4.50 


-0.17 


-2.4 


2.37 


giant planet, Lup IV 


3 




-10.1 


5.6 


1000 


8.00 


-2.18 


-2.6 


0.39 


giant planet, Lup III 


4 




-9.7 


< 1.6 


600 


4.50 


-1.05 


-2.2 


0.30 


grain growth 6 , Lup III 


5 




-10.6 


< 3.4 




2.20 


-0.93 


-1.4 


0.86 


grain growth 6 , Lup III 


6 




-9.5 


< 1.8 


140 


5.80 


-1.42 


-2.6 


0.43 


grain growth, Lup III 


7 




-9.1 


1.1 


80 


8.00 


-1.55 


-2.9 


1.27 


grain growth, Lup III 


8 




-11 


1.6 


76 


5.8 


-0.86 


-2.8 


0.17 


grain growth 6 , Lup III 


9 




< -11 


< 3.4 


560 


4.50 


-0.31 


-1.9 


0.22 


photo, disk, Lup III 


10 




-9.3 


5.1 




2.20 


-1.19 


-1.5 


0.27 


grain growth 6 , Lup III 


11 




-9.6 


2.7 


230 •■• a 


4.50 


-1.05 


-2.3 


0.96 


circurnbinary /gr-grow 6 , Lup III 


12 




-11 


3.8 




5.80 


-0.28 


-2.1 


0.11 


grain growth 6 , Lup III 


13 




< -11 


< 0.8 


234,390" 


8.00 


-1.67 


-3.7 


0.48 


circumbinary/debris, Sep 


14 




< -11 


< 1 


20 


5.80 


-0.42 


-2.7 


0.14 


circurnbinary/ photo, disk, CrA 


15 




-8.6 


< 2 




8.00 


0.76 


-2.4 


0.46 


giant planet, CrA 


16 




< -11 


< 2 


75 


8.00 


-1.74 


-3.2 


0.07 


debris disk, CrA 


17 




-11 


< 2 


75 


5.80 


-1.06 


-2.3 


0.06 


grain growth, CrA 



a Objects # 11 and 13 are triple systems. The distances given correspond to the closest and widest components. For the case of # 11, 
the tightest components can not be resolved. 

b Uncertain classification due to relatively weak evidence for accretion (target # 8 and 12) or SEDs similar to those of classical TTauri 
stars of spectral type M (target # 4, 5, 10 and 11). 
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Table 7: Fractional distribution of transition disk organized by clouds 



Disk candidates 


Lupus III, IV 


Cra 


Opht 


debris 


- 


1 (25%) 


4 (15%) 


photoevaporated 


2 (17%) 


1 (25%) 


5 (19%) 


grain growth 


8 (66%) 


1 (25%) 


13 (50%) 


hosting giant planets 


2 (17%) 


1 (25%) 


4 (15%) 


circumbinary 


1 (8%) 


1 (25%) 


4 (15%) 


Total 


12 


4 


26 


Age [Myr] 


~ 1.5-4 a 


~l fe 


~ 2 C 



t From Paper I 



References: a Comeron (2008); b Wilking et al. (2005); c Sicilia-Aguilar et al. (2008) 



